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We solve for the exact atom-held eigenstates of a single atom in a three dimensional 
spherical cavity, by mapping the problem onto the anisotropic Kondo model. The spectrum 
has a rich bound state structure in comparison with models where the rotating wave 
approximation is made. It is shown how to obtain the Jaynes-Cummings model states in 
the limit of weak coupling. Non-perturbative Lamb shifts and decay rates are computed. 
The massive Kondo model is introduced to model light localization in the form of photon- 
atom bound states in photonic crystals. 
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1. Introduction 


The system of two-level atoms coupled to radiation is a fundamental problem of 
physics, and much work has been devoted to its study, especially in the context of quantum 
optics|L|. For the most part, a satisfactory understanding of the problem is obtained after 
making some approximations, most importantly the rotating-wave and slowly varying en¬ 
velope approximations. Though these approximations are well-justified near resonance, it 
is important to have an understanding of what is missed in making these approximations, 
and to ascertain whether there are novel effects of any significance. 

In this paper we present an exact solution to the problem of a single two-level atom in 
a spherical cavity, without making the above approximations. This is a considerably more 
complicated problem because of the vacuum fluctuations. The solution is accomplished 
by mapping the optical problem onto the Kondo problem of massless fermions interacting 
with a magnetic impurity, as was done for a one-dimensional fiber geometry in |2j|. 

The novel features of our solution, which are unanticipated based on previous studies, 
have to do with the spectrum of eigenstates of the atom-field interaction. For a single 
cavity mode in the rotating-wave approximation, the atom-field eigenstates are known 
from the solution of the Jaynes-Cummings model. The parameters of the latter model are 
the two-level splitting cu 0 , and the coupling a of the atom to the radiation. The latter 
determines the decay rate Tdecay = o?L where L is the volume. The eigenstates of the 
Jaynes-Cummings model correspond to doublets of atomically dressed superpositions of 
N and IV — 1 photon states. Our solution is given below in terms of a scattering theory 
description for a spectrum of massless particles, and we describe how information about the 
atom-field eigenstates is encoded in the scattering matrices. As a result, we find that the 
spectrum of states in our model is much richer, in a way that depends on the dimension¬ 
less ratio g — T^ ec&y / ttojo, which is a quantum coupling. Roughly speaking, we find that 
the one-photon (IV = 1) eigenstates are 1/g in number, the lightest becoming identified 
with the Jaynes-Cummings one at weak coupling (small g) and near resonance, and the 
higher ones arising as bound states, and as solitons. Interestingly, at the special strong 
coupling point g — 1/2, the fundamental Jaynes-Cummings eigenstate becomes unbound 
and disappears entirely from the spectrum, leaving only solitonic states. In addition to 
these spectral features we are able to compute non-perturbative expressions for the Lamb- 
shifted energy splittings and decay widths. 

The above spectrum has a sine-Gordon-like character. Sine-Gordon theory in not 
unfamiliar in this context, as it is well-known that a resonant dielectric medium of two-level 
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atoms in a semi-classical limit is well-described by the classical sine-Gordon theory. This 
is the phenomena of self-induced transparency0. Our work can be interpreted as showing 
that the propagating spectrum in the theory of self-induced transparency has an important 
significance for a single atom in a cavity interacting with fully quantized ratiation. Our 
theory thus provides a bridge between the theory of self-induced transparency and the 
Jaynes-Cummings model. 

In the final sections of the paper we apply our techniques to the problem of an atomic 
impurity in a medium with a photonic bandgap, with applications to photonic crystals in 
mind0. Photonic crystals, envisaged in early works of Yablonovich and Jolm[^J |J, are 
periodic dielectric structures exhibiting gaps in the allowed energies of photon propagation, 
in close analogy to electronic band structure. Important effects are expected when these 
materials are doped with atomic impurities. Specifically, the spontaneous emission of an 
atom is severely inhibited if the energy splitting of levels in the atom coincides with an 
energy that is forbidden to propagate in the material. Under these circumstances, the light 
can form a bound state with the atom, this photon-atom bound state being the optical 
analog of an electron-impurity-level bound state in the gap of a semiconductor]]]]. In this 
paper we introduce a toy model where these effects can be studied exactly by adding a 
mass term to the Kondo model in a way the preserves the integrability of the latter. Here, 
the origin of the photonic bandgap is the gap between polariton branches. We solve this 
massive Kondo model at a special point where it is equivalent to a free fermion theory, and 
indeed find a photon-atom bound state with the expected properties. We can also compute 
exactly the binding energy of this state. A different integrable model of this phenomenon 
was studied in § wherein the rotating wave approximation is made. 

We present our results in the following way. In section 2 we use spherical symmetry 
to reduce the problem from three spatial dimensions to one, and in section 3 describe 
the resulting theory as a boundary quantum field theory. For the resulting theory to be 
exactly solvable, one needs to make an approximation that favors photons in the vacinity 
of the resonance; such an approximation is not necessary in the one-dimensional case. In 
section 4 we map the problem onto the anisotropic Kondo model, describe the scattering 
spectrum in infinite volume, and present exact reflection S-matrices expressed in terms of 
the Lamb-shifted resonant energies and decay widths. In section 5 we describe the meaning 
of topological charge in the optical context. In section 6 we describe how information about 
the atom-field eigenstates is encoded in the reflection S-matrices by computing the one- 
particle finite volume spectrum and comparing this to the Jaynes-Cummings spectrum. 
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Here we show how for each particle of the infinite volume spectrum one can associate 
polariton-like states with optical phonon behavior at large and small energies. In section 
7 the massive Kondo model is introduced to model atomic impurities in photonic crystals. 
The latter model is solved at the special point g — 1/2 in section 8, and this solution is 
used in section 9 to study the localization of light. 


2. Spherical Dimensional Reduction 

We consider a two-level ‘atom’, with spherically symmetric states |±) which are eigen¬ 
states of the unperturbed atomic hamiltonian // ( 'j l,om : 

«T m l±> = ±y l±>, (2.1) 

such that the energy splitting is uoq. We work in spherical coordinates r, 9 , 6. and place the 
atom at r = 0. We will focus on electric-dipole transitions, thus we further assume |—) is 
an angular momentum eigenstate (/, m = 0), whereas |+) is an (l + 1, m = 0) state. The 
wave-functions of these states are then 

<r1+> = f+(r)Y l+ i i0 (M), <fH = (2.2) 

In coupling the atom to radiation, one uses the electric-nrultipole expansion as de¬ 
scribed in e.g. j9]. Expanding electric and magnetic fields in frequency components, 

E(f, t) = Jdw e~ iujt E u {r) + e^E^, (2.3) 

and similarly for B. and keeping only the electric-dipole term in the multipole expansion, 
one has 

= 2 uo 3/2 a(u)ji(ur)X 10 (9,(/)) 

(2.4) 

Eu(r) = 2 iu: 1/2 a{u) V x j\{ur)X 1Q {p, 4>) 

where X/ m = LY[ rn / \Jl[l + 1). and ji is a spherical Bessel function. Using orthogonality 
relations for Xi m and ji, the free field hamiltonian is 

tffieid = J d 3 r^E 2 + B 2 ^ — duo uoa\uo)a{uo). (2.5) 

In the quantum theory one has the commutation relation 

[a(a;), a){u')\ — 5{uo — uo'). (2.6) 
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The atom is coupled to the radiation with the usual dipole interaction H int — —d ■ E 
where d is the electric dipole operator of the atomic electron, d — er. Letting r — rr, one 
can show 

r ■ E u (f) = — a/8c J a{uj) 


10 - 


(2.7) 


We assume that the wavelength of relevant photons is large compared to atomic di¬ 


mensions, such that the r dependence in (2.7) can be pulled out of the matrix element: 

-(+\d-E u (r)\—) m lim a/ScJ a (w)(+|erY 10 |-). (2.8) 


r —>0 


The matrix element (+|erYi 0 |—) can be related to the usual reduced dipole moment 
d defined as (+|d m |—) = dS m o, where d m is the m = ±1,0 spherical component of the 


vector d. Using r m — r\j one finds 


(+\erY 10 \-) = -^d, 



where 


d = e- 


(/ ± 1 ) 


drr 3 /*(r)/_(r). 


(2-9) 


( 2 . 10 ) 


yj(21 + 3)(2 1 + 1) J o 

Finally, using linx^o ji(x)/x — 1/3, and the reality of Yjq. the complete hamiltonian 


is 


rr _ rrfield . jratom . tt 

H — + Uq + Hi nt 


where 


and 


H t 


atom ^0 


0 


= f 5 3, 


Hint = y — { j duju 3 / 2 (a(uo)e- iujt + a\uu)e iujt )} (d S + + d* S ~) . 
Here, S 3 , S' 1 * 1 are Pauli matrix operators: 

[S 3 .S ± ] = ±2S ± . [S + , S~] = S 3 . 


( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 


(2.14) 


3. Boundary Field Theory Description 

In this section we will give a boundary field theory description of the model of the 
last section. Boundary field theories live on the half-line in space, x > 0, with non-trivial 
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interactions at the boundary. The half-line occurs automatically in the spherical reduction 
since r > 0, and the boundary interaction is at r = 0. 

The free hamiltonian Illj cld is equivalent to a free boson on the half-line: 

H ° eu = r dr \ + (9r ' # ’ )2 ) ’ 


with the supplemental Neumann boundary condition d r (j){r — 0. t) = 0. The mode expan¬ 
sion of (f) can be written as 


<t>(r, t ) = i 



dk 1 




(3.2) 


where k ■ x — \k\t — kr. The Neumann boundary condition enforces ci(k) — a(—k ), so that 
at r = 0, cf)(r, t ) can be expressed as an integral over k > 0 modes. 

In order to make the problem solvable, as we will see, we make an approximation that 
favors photons in the vacinity of u> o, so that cn 3 / 2 in (|2.13|) is replaced by Wqa/uJ. Then, 


duj to 3 / 2 (ci(uj)e luJt + (cu)e' lU}t ) « y/iTuJodt(j)(0,t). 


(3.3) 


If d is complex we absorb the phases of d. d* into the definition of S± without changing 
the commutation relations (|2.14). Finally, we obtain the hamiltonian! 


H = H™ d + + jdt<K 0. t) (S+ + S~), 


(3.4) 


where 



(3.5) 


To lowest order in perturbation theory, the hamiltonian ( |3.4|) gives the well-known result 
rdecay = | d 2 u>Q. Our solution below presents all-order (d 2 corrections to this result based 


Though we need to make the approximation ( |3.3|) in order to solve the model, for¬ 
tunately in a one-dimensional fiber geometry, one finds u; 1 2 rather than cn 3 / 2 in ( |2.13| ), 


2 This hamiltonian is closely related to the so-called Lee model and its generalizations] 10], 
except that in our model the fermion fields are fixed at one spacial location. The latter two 


papers in ||10| study the renormalization problems associated with what amounts to the rotating- 
wave approximation. 
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so that the model is exactly solvable without any further approximations [pi]. All of the 
subsequent results of this paper hold in a fiber geometry with 


/3fiber 




(3.6) 


where A e ff is the effective cross-sectional area of the fiber. 

The parameter /3 2 /8ir is a dimension-less coupling constant, where the strong coupling 
regime is large /3 2 . Using the lowest order perturbative result for T^ecay, this coupling can 
also be expressed as 


9 = 




i r 


decay 


(9 « 1). 


(3.7) 


8tT 7T Ido 

In order to get some idea of the size of the coupling for real atoms, we can consider a 
hydrogen-like atom with the nuclear charge Ze. The largest value of (3 occurs for the Is to 
2 p transition, where g — Z 2 e 6 2 11 3 _9 /7r ~ Z 2 10 -8 . Here, g goes as e 6 since the Coulomb 
interaction determines both ujq and d. One can perhaps hope for larger values of g in 
artificial atoms, e.g. quantum dots, wherein d and ujq are fixed by distinct physics. 


4. Mapping to the Kondo Model 

4-1. The Exact Scattering Spectrum 

In |Q], a single atom in a fiber, which is described by a hamiltonian of the form (53) 
except that the theory lives on the full line —oo < x < oo and j3 is different, was mapped 
onto the Kondo model. Here, the result is even simpler since one doesn’t have to fold the 
system. Namely, H is related by a unitary transformation U to the bosonized form of the 
anisotropic Kondo hamiltonian Hk'- 


H = U j H K U , 


(4.1) 


where 


and 


h k = H$ eld + y ^ s + e wm/2 + s - e -n3m/^ 5 


u 


J_ i/3Ss<M0)/4 

V2 


(S 3 + S+ + S-). 


(4.2) 


(4.3) 
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Since H and II k are related by a unitary transformation, the quantum mechanics can 
be formulated in either the ‘optical picture’ or the ‘Kondo picture’. Matrix elements of 
operators O in the two pictures are simply related: 


I ^ IV 7 ) optical = K^'\OkW)K, 


(4.4) 


where | i/j)k — U\i/j), and Ok — UOU^. We remark that the trivial atomic hamiltonian 
Hq to™ j n the optical picture becomes the interaction in the Kondo picture. 


Both models ( |3 . 4| )( |4.2| ) have ultra-violet divergences, and furthermore the quantum 
operator U is in need of regularization. We will take the point of view that the equation 
( |4 . 1| ) defines the regularization of H once we have regularized Hk and U. The regulariza¬ 
tion of Hk is as in the sine-Gordon model and amounts to properly normal-ordering the 
exponential operators; this leads to the anomalous scaling dimension —g for the parameter 
ujq. It can be checked that this is consistent with the renormalization performed directly 


in the optical hamiltonian H as was done in |JTT[| , where it was shown that to lowest order 
the beta function reads gd^cuo — —gwo- 

The physical parameters of the optical problem are g, which is governed by the 
strength of the dipole coupling, and the energy scale set by the two-level splitting. The 
parameter wo in ( |4.2| ) is a bare, unphysical parameter; physical energy scales are a function 
of Wo, g and an ultraviolet cut-off /i. In the scattering theory solution presented below, 
physical energy scales will be set by the parameter wg. This scale is a known function 
which we will not need, but for completeness include from [El: 


UJb — 


Toot (W»- 

7T \2 — 2g 


1- 2 g 

2- 2 g 


2-3 g 
2-2 g 




V 2 


(4.5) 


In the Kondo model the physical parameters are the dimension-less anisotropy pa¬ 
rameter g and the ‘Kondo-temperature ’ t k , which is related to u>b as follows: 


Tk = tan 


-Kg 


ujb- 


(4.6) 


, 2-2 g, 

Tk is defined such that the atomic impurity contribution to the partition function behaves 


as 


•^atom — 2 cosh (: Tk/T ) + O(g), 


(4.7) 
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thus Tk is convenient for describing thermodynamical properties, as was done in |2|. For 
small g, Tk ~ u>o/2; this is consistent with the fact that when g — 0, Z A tom should 
represent the partition function for a two-level system with energies ±u>o/2. Henceforth, 
we will express all energies in terms of ujb- 

Most importantly, the spectrum of H and Hk are identical. Thus we can infer what 
are the quantum states that exactly diagonalize the atom-field interaction from knowledge 
of the exact spectrum in the Kondo modelS. This spectrum consists of a rich spectrum of 
massless particles which have a sine-Gordon-like character. The origin of the sine-Gordon 
spectrum is due to the fact that a bulk sine-Gordon interaction is compatible, as far as 
integrability goes, with the boundary interaction. (See below.) This is analogous to the 
treatment of the boundary sine-Gordon theory in |T5[, the difference being in the reflection 
S-matrices. Namely, the spectrum consists of a soliton and anti-soliton, and [1 jg — 2] 
breathers which can be viewed as soliton-anti-soliton bound states. The meaning of this 
spectrum will be elaborated upon in section 5. The point g — 1/2, the so-called Toulouse 
point of Kondo physics, can be formulated as a free fermion theory using bosonizatiou. 
(See below.) In the repulsive regime g > 1/2, there are only solitons and anti-solitons. For 
g > 1, the interaction is an irrelevant operator and the theory breaks down. In the bulk 
these particles propagate with a massless dispersion relation: 

E a = P a — Ha e ° f° r right — movers 
E a = — P a = g a e~ e for left — movers (4.8) 


where a G {1, 2,.. < (1 — g)/g, s, s} is an index running over the breathers and solitons, 
and 6 is a rapidity variable. The parameters g a are related as follows: 


Hn 


2 g sin 


( n-ng \ 
V2-2 g) ’ 


Hs Hs — Hi 


(4.9) 


where g is an arbitrary, unphysical energy scale which can be set to 1. 

The interactions of these particles with the atom are encoded in the exact reflection 
S-matrices for these particles at the boundary [|I2|| [|T3| . Let R b a (9 ) denote the S-matrix for 


6 We refer the reader to [12] |13] for a description of the Kondo model that is the most useful 
in our context. An exact Bethe-ansatz solution of the Kondo model was found in [14]. References 
to other earlier papers can be found in [13]. 









particle of type a to scatter off the boundary and be reflected into a particle of type b. 
Defining 0 B as: 

uj B = fie eB , (4.10) 


the explicit expressions for the reflection S-matrices for the left-movers are 


R+m 

Rim 


R-(9) = tanh ^— ~2~~ 
RZ(9) = 0 


Rn( 0 ) = 


tanh 


nrgn 


tanh 


( e+e B 


+ 


vngn 

4(1^) 



(4.11) 


Here + = s and — = s. 

It will be useful to express the reflection S-matrices in terms of resonances u5o, n? widths 
r n , and physical energies E: 


R+(E) 

Rn(E) 


+ E 2 

So,s - E 2 + ir s E 
Wo, n " E 2 -iV n E 
ul n -E* + iT n E' 


(4.12) 


where 


^o,s = wb, r s = 2 uj b 


^0,71 


2 sin 


( n-ng \ 
V2-2 g) 


wb, 


4 sin 2 


/ ring \ 
V2-2 g) 


UJ B - 


(4.13) 


Near the resonances E rs u)o,n> the reflection S-matrices have a Lorentzian signature: 


Rl 


(r„/ 2 ) : 


(Do.n - E ) 2 + (r„/2) : 


For the solitons, one has 


r: 


(T,/ 2 )- 


(So,. - E) 2 + (r,/2) : 


(4.14) 


(4.15) 


As we will show in section 5, the quantities cDo. n . F n determine the energy spectrum of 
eigenstates in finite volume. 

It is well known in the quantum sine-Gordon theory literature that the n — 1 breather 
is the particle corresponding to the scalar held (j) itself [|T6f . To make this more explicit, in 
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[|] electric field correlation functions were computed in the single 1-breather approximation 
which is remarkably good for values of g up to ~ 1/5. Because of ( |4 .1 4| ) , the vacuum power 
spectrum is Lorentzian where Dq.i and Fi represent the Lamb-shifted two-level atomic 
energy splitting and decay rate. These agree with lowest order perturbative computations: 


w 0) i = uj 0 + O(g) 

Ti = gnuo + 0 (g 2 ). 


(4.16) 


5. Topological Charge 

In the usual bulk sine-Gordon theory on the full line, the soliton and anti-soliton carry 
topological charge Q = ±1, where 

Q = if dx d x (f) — (4>(oo) — (f>(—oo)) . (5.1) 

In our boundary version of the problem, solitons are still characterized by a charge Q — 
/3(4>(r — oo) — 4>{r = 0))/27r. The value of (f> at r — 0 isn’t fixed, thus topological charge 
conservation can be violated. For simplicity, set <f)(r — cx)) = 0. An incoming soliton 
(' Q — 1) with (f){r — 0) = —2tc/(3 can be reflected at r — 0 into an outgoing anti-soliton 
with <j)(r — 0) — 27 r//3, violating charge conservation by 2 units. In this process, the field 
at the boundary interpolates between —27 r//3 and 2 tt// 3. The fact that /?z = 0 means that 
only processes which violate topological charge conservation by two units are allowed in 
the soliton sector. 

We can relate the topological charge of solitons to simple properties of the electric 
field E oc dt(j). Imagining the soliton as a pulse localized in space, then an incoming soliton 
at time —t and an outgoing anti-soliton at time t. where t — 0 corresponds to complete 
absorption, have pulse profiles in space simply related by (f) —> —0. Thus, the phase of 
the electric field of the outgoing anti-soliton is shifted by 7r in comparison to the incoming 
anti-soliton. 

Consider a soliton scattering process where the atom is initially in its ground state. 
When the soliton is absorbed, the atom is precisely in its excited state; it then emits an 
anti-soliton and returns to its ground state in the far future. This can be seen by examining 
the original optical states |±) in the Kondo picture. Defining |±,(/>(0))x = U |±), one has 

\±,<j>(0)) K = (e w(0)/4 |+) ± e -W(o)/4|_^_ (5.2) 
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Note that 


|“» 0 (°) + -p)K = i\+,<t>{fy) K , \-,(/){0) + —) K =-\~,(/)(0))k- ( 5 . 3 ) 

Thus, a soliton incident on the atom in its ground state in the far past corresponds to 
|—, —27 t//3), and when the soliton is completely absorbed 0(0) = 0 and the state has 
evolved to i|+, 0), i.e. is in its excited state. In the far future, 0(0) = 2n/f3, and the 
state is — | —, 27 x/j3). The atom has necessarily returned to its ground state due to energy 
conservation. Note however that the atomic state has changed phase by 7r. 

The fact that soliton absorption excites that atom to precisely its excited state in¬ 
dicates that the solitons are the most fundamental excitations. As we will see below, at 
strong coupling only the solitons remain in the spectrum. 


6. Optical Phonon Spectrum and the Jaynes-Cummings Model Limit 

Though the scattering theory described above provides an exact solution of the model, 
its precise meaning is not entirely clear from what we have done so far. The original 
Hilbert space, TL. consisting of the two-level space tensored with the photon Hilbert space, 
has been replaced with an infrared description 7i pa , rticle- The latter Hilbert space contains 
a spectrum of particles with continuous energies, and does not have a two-level structure. 
The distinguishing feature of this spectrum of particles is that they provide a basis that 
diagonalizes the atom-field interaction; the existence of the two-level atomic Hilbert space 
is encoded in the reflection S-matrices R. A similar treatment of the simpler case of 


harmonic-oscillator atomic impurities was given in [17]]. Irrespective of these remarks, the 
free field hamiltonian I/),' 01 ' 1 is a linear theory in vacuum, so an intrinsic understanding of 
the non-linear sine-Gordon-like spectrum seems called for. 

It is also desirable to describe eigenstates of the atom-field interaction in a fashion 
that is more understandable in terms of the original Hilbert space 7i. This is a difficult 
problem in general, since it requires computing the inner-product of states in 7i with 
states in 7T r)ar t.icic- Exact eigenstates can be described in the space 7 i for the much simpler 
Jaynes-Cummings (JC) model. In this section we describe how information about the 
atom-field eigenstates is encoded in the scattering theory description, and indeed we show 
how to obtain the known JC eigen-energies in the appropriate limit. This exercise leads to 
the understanding that for each particle of the spectrum of section 3 there is an associated 
optical phonon, in a sense that will be explained. 
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6.1. Jaynes-Cummings Model 


The Jaynes-Cummings model [|I 8 | is an exactly solvable model of the atom-field inter¬ 
action which arises under various approximations to the hamiltonian (|3.4| ). Let us suppose 
that the one-dimensional volume of the cavity is L. The first approximation is to only 
consider a single discrete photon mode of the cavity of frequency co c near ta 0 . Namely, 
the free photon Hilbert space is built from the single mode creation-annihilation operators 
a = a{uo c ) and ad = ad(ca c ) satisfying [a, a'*'] = 1. Next, we make the rotating-wave ap¬ 
proximation and drop the terms aS~ and a* S + in (|3.4|) . Finally, we drop the exponentials 
exp(ca c — Wo) in the interaction. Recalling that in finite volume f dk/'/Jn —> one 

obtains the hamiltonian: 


H — oo c ad a + —~ S 3 + a ( aS + + ad S ) , 

Zj 


where 


a 


nguic 


( 6 . 1 ) 


( 6 . 2 ) 


LUQ 


The volume L drops out of many physical quantities; for instance T^ecay = ct 2 L| Wc= 

(3 2 u) 0 / 8 . In the limiting case of a micro-cavity of size tuned to the resonant frequency 
ui c — to 0 with uiq / 7 r = 1 /L, one has a m i cro = y/gnw o- 

The above hamiltonian is easily diagonalized exactly. Let \N, ± > denote a basis of 
states with N photons and ± the two-level states of the atom: 


a'a|JV, ±) = N\N, ±), S 3 |JV. ±) = ±|JV, ±). 


(6.3) 


The hamiltonian has a 2 x 2 block diagonal form where the blocks act on the pairs of states 
|IV — 1, +), |IV, —). Each block is easily diagonalized leading to an energy spectrum £^, 


N = 1 , 2 ,...: 

= w c(IV - 1/2) ± ^(Wc-Wo) 2 +4ct 2 IV, 
where the eigenstates are given by 

= cos ©|IV — 1, +) + sin ©|IV, —) 

= — sin ©|IV — 1, +) + cos ©|IV, —), 


(6.4) 


(6.5) 


where 

f) _ ■>/( w c — w o) 2 + 4a 2 IV + ( uj c — cuo) 
y/{u c - wo) 2 + 4a 2 IV - (u c - u 0 ) 
This spectrum is observable experimentally ||19f . 


( 6 . 6 ) 
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6.2. Optical Phonon Spectrum 

Let us describe how to obtain the JC spectrum from our scattering theory under 
certain limits. In order to compare with the single cavity mode approximation, we must 
place the atom inside a finite volume, letting the spherical cavity have a diameter L. There 
is now a quantization condition on allowed momenta which depends on the reflection S- 
matrices. Consider a wave-function for a particle from the spectrum of section 3 as it 
traverses the radius of the cavity, is reflected and returns to the edge. This wave-function 
picks up a phase which must equal unity: 

e iPL R{E ) = 1. (6.7) 

Taking the logarithm, and using P — — E , one obtains 

E+UogR(E) = ^, ( 6 . 8 ) 

where m is some integer. The quantity on the RHS determines the mode of the cavity: 
oo c = 2tc m/L. 

Let us now make two approximations. As in the JC model we assume the cavity 
is tuned near resonance uj c ~ <2'o. n for the quantization of the n-th breather with R = 
R™ in (|OD. Secondly, at weak coupling, for the low n breathers, it is safe to suppose 
T n -C |cDo,n — E |, since T n ~ n 2 n 2 g 2 ujB and cDo ,n ~ nngujB- Under these approximations, 
logi?™ ~ — zT n /((Lo,n — E). Inserting this into (|678|), one obtains a quadratic equation for 
E , the solution for the n-th breather being: 

1,„ W = —— ± l (K - Co ,„) 2 + 4T„/L) 1/2 . (6.9) 

Multiparticle (N > 1) eigenstates are characterized by more complicated quantization 
conditions involving the bulk S-matrix ||15|| . Since the ground state energy of the JC model 
is — u>o/2, whereas that of the scattering description is set to zero, one should subtract 
ujo/2 from the above energies in order to compare. Using the small g limits, 

E N=i,n( u c) - uq/2 S3 ^ + UJ C ± i ((w c - nuj 0 ) 2 + 4n 2 gnuJo/L) 1/2 . (6.10) 

From ( |6.2| ), with uj c ~ loo, one sees that for the n — 1 breather, one-particle energies 
correspond precisely to the one-photon JC energies: 

En=i,ti=i — w o /2 = ( 6 - 11 ) 
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At oo c — ooq, the N = 1 JC eigenstates are symmetrically split about ooq/2. Let us 
define then a scaled energy splitting E' such that: 

~ u 0 /2 = ±E , = ±( g / 2 )i/2 m ( 6 . 12 ) 

U >0 

We can then compare this with the ‘exact’ prediction (|6.8|). Setting E' — (E — cDo,i)/^o,i 
in ( |6.g| ) where R is the reflection S-matrix for the first breather R \, and oo c — cDq.i = 27r/L, 
one finds that E' is a solution of the equation 




E' 1 + 2 E' + 2 i sm(ng/ (2 - 2 g)){E' + 1) 
E' 2 + 2 E' — 2 i sin(7T(7/(2 — 2g))(E' + 1) 


= 0. 


(6.13) 


The differences between the JC model and the result (|6.13|) are displayed in figure 1. At- 
small g they of course agree very well. The strong departure from the JC model at g — 1/2 
is due to the fact that this is the threshold for the disappearance of this particle. 

The higher n-breather N — 1 particle states are not predicted by the JC model. We 
believe this is because higher resonances E ~ uojq are eliminated in the rotating wave 
approximation. Based on the above observations we can formulate the following appealing 
picture of the spectrum. At large volumes, the energies behave as 


-®iV=l,n ~ u 


K=i. 


Cl 


^0, 


E 


N=l,n 

E n =i, 


W0,r 


00 r 


LO c ^ 0 , 77-5 

LO c ^0,n* 


(6.14) 


Thus the E + branch is optical phonon-like at small uo c and photon-like at large oo c , whereas 
for E~ this is reversed. For n — 1, these states are the precursor to the well-known polari- 
ton branches which occur in the context of a dielectric medium of atoms. (See for instance 
and the next section.) Here, there is no medium, however these ‘vacuum-polaritons’ 


are characterized by the regimes where they are phonon-like, so we will simply refer to 
them as optical phonons. The n — 1 breather optical phonon is the basic excitation of the 
JC model. Since the n-th breather leads to an optical phonon at energy approximately 
nooo for small g, the latter can be understood as a bound state of n optical phonons of 
energy u 0 . 

We thus propose that our model contains a rich spectrum of optical phonons of energy 
cDo, n , cDo,s: in correspondence with the breather/soliton spectrum in section 3. The n-th 
optical phonon is a bound state of n fundamental n — 1 phonons. As g increases the 
higher optical phonons become unbound and disappear from the spectrum one by one. 
When g — l/(n + 1), from ( |4.13| ) one sees that Qo,n — 2cDo ,s- Thus as g increases and 
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Figure 1. Comparison of the energy splitting of the N — 1 photon Jaynes- 
Cummings eigenstates and the predictions of our model. The y axis is E , 
which, for the JC model is defined in ( |6.12| ), whereas for our model it is a 
solution to Q6.13 ). 


reaches the point l/(n + 1), the n-th optical phonon unbinds into two solitonic phonons. 
This also implies that all of the n-th breather optical phonons can also be viewed as a 
bound state of two solitonic optical phonons. Finally, when g — 1/2, the n — 1 optical 
phonon, which is the fundamental excitation of the JC model, also disappears from the 
spectrum leaving only the solitonic phonons for all g > 1/2. 


7. Atomic Impurity in a Photonic Crystal 

In this section we place the atom in a dielectric medium. One can generally model 
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such a medium with a dielectric constant 


£ ( UJ ) — £ O0 + ( £ 0O — £()) n- T 2 5 (7-1) 

L0 Z — Wy 

giving the dispersion relation e(lo) = k 2 )co 2 . (See for instance 0.) For simplicity, let us 
set £00 = 1. The dispersion relation has two polariton branches co±(k) as shown in figure 
2 . Since co+(k = 0) = and cu-(k = 00 ) = lot, there is a gap between the two 

branches E gap — (^/Eq — 1 )lot- Though this is not the manner in which photon bandgaps 
are thought to arise in e.g. 0, we believe it serves as a good model for the physics we are 
trying to study. 



Figu re 2. The two polariton branches of the dispersion relation following from 
(7.1), with £0 — 3. Both axes are in units of lot■ The upper branch intersects 
the y-axis at y/Eg. 
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In a certain limit, the dispersion relation is approximately relativistic. Namely, let us 
treat the lower branch as constant u>-(k ) = ut, and shift the energy of the upper branch 
u;+ = cut + c*/. Then, if £q is large, w' is large compared to ut and u/ 2 — k 2 ~ m 2 ap with 

m gap = \/£o - 1 w T . (7.2) 


For £ 0 1, m g a p ~ E gap . Alternatively, one can view this limit as taking ut —> 0, and 
£q —> oo, keeping m gap = yGo wt fixed. In this limit the lower branch disappears and the 
upper branch is a massive relativistic dispersion relation. 

One can model this gap in the medium while preserving the integrability by adding a 
sine-Gordon interaction in the bulk: 


H, 


field 


H, 


field 


m 


-jp J drcos(/?0). 


(7.3) 


Choosing m as in (|7.2|) , the above interaction gives the photon the proper dispersion 
relation. Mapping onto the Kondo model as in ( |4.1|) , one now obtains 


H 


i m ) _ 

K ~ 


dr 


l 2 (d t 4>) 2 + ^( d r <j> ) 2 - cos (P<f>) 


+ 


— ( S + e' w(0)/2 + S~e~ i ^ ( - 0)/2 

2 V 


(7.4) 


The desirable feature of the theory (HD is that the bulk and boundary interactions are 
compatible as far as integrability is concerned. Namely, the exact bulk sine-Gordon spec¬ 
trum also diagonalizes the boundary interaction with the atom. The spectrum of the model 
is the same as in section 3, except that the particles are now massive with a mass scale 
determined by m gap . The dispersion relation is now 


E a = M a cosh 0. P a = M a sinh 6 . (7.5) 

satisfying E 2 — P 2 — M 2 . The masses M a are related as in ( ff~9l) , 

M n = 2 M s sin , (7-6) 

but now the soliton mass M s is a physical parameter depending on m gap . Since the first 
breather is identified with the photon, we will set 


Mi = m gap . (7.7) 

In the literature on the Kondo model, a bulk mass term as in ( |7.4|) is never considered 
since one deals with massless fermions near the Fermi surface. We believe the model 
( [7.3| ) is integrable since the bulk massive spectrum is compatible with the spectrum that 
diagonalizes the boundary interaction. The boundary reflection S-matrices for the model 
( [7.4| ) are unknown. For simplicity, in this paper we will limit ourselves to the model at 
g — 1/2 and derive the S-matrices in the next section. We intend to present the general 
case elsewhere. 
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8. Free Fermion Point: Massive Case 


Remarkably, at the point g — 1/2 in the strong coupling regime, the theory is equiv¬ 
alent to a free fermion theory. Though this is well-known in the Kondo physics literature, 
this is a novel phenomenon in the optical context since at this value of the coupling the 
perturbatively dressed photon, in the sense of the Jaynes-Cummings model, disappears en¬ 
tirely from the spectrum. The theory at this coupling will be considered in greater detail 
in ]2T[] . Here, we consider the bulk massive case, which has not been studied before, with 
applications to photonic crystals in mind. 


8.1. Physical Parameters 

We first define the physical parameters. When g — 1/2, Tk diverges however cub 
remains finitcS. For g — 1/2 we define the physical parameters as 


ujq = 2 ( Uo, s = 2 ujb, 


decay 


2(Uq. 


( 8 . 1 ) 


Note that this precisely corresponds to the 1-breather expressions cDq,i and Id evaluated 
at g = 1/2. The parameter To represents the energy splitting of the atom; note that since 
this is twice cDo, s , the soliton is in resonance with half of the energy splitting. 

From (7.0) one sees that the breather mass M\ is exactly twice the soliton, or fermion, 
mass M s , which reflects the fact that the lowest breather is a soliton/anti-soliton bound 
state and just becomes unbound at g — 1/2. In this section, we will express everything in 
terms of the fermion mass M = M s — m gap /2. 


8.2. Reflection S-matrices 

We now derive the reflection S-matrices in a treatment that is similar to what was 
done for the boundary sine-Gordon theory in |[T| . 


In the massless limit, the scalar field decomposes into left and right-moving compo¬ 
nents (f) — tp{z + ) + where z± = t ± r. The Neumann boundary condition d r <f> — 0 

then reads <9 z +y?(z + ) = d z -Tp{z~) at r = 0. This implies 


p{z ) = P( z ) 


a 


\f\n r ’ 


( 8 . 2 ) 


where a is some constant. 


One finds cub oc cUq, which is consistent with the anomalous ydnass dimension of coo- 
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The fermion fields, with topological charge ±1, are the following: 


= e ±iVI ” ip ± = e TiV ^ tp . ( 8 . 3 ) 

In terms of the fermions, the Neumann boundary condition reads 

— e Tlcr (r = 0) (8.4) 

which breaks the topological charge symmetry. The action which enforces the Neumann 
boundary condition is 

/•OO 

Sfree = / drdt [iip-(d t - d r )ip + + #_(< % - d r )ip + - + ip-ip + )] 


(8.5) 


— ij dt (e l<J '0 + '0 + + e . 

Using (|8.2|) , the interaction can be written as 

5 in t = ~ J dt [(ip+a- + ip_a + )S + + 5_(a+^_ + a_V’+)] , (8.6) 

where a± — e^ ia l 2 . The parameter A, with units of yhnass, will be related to below. 

The a± should be regarded as fermionic operators: 


a±ip — —ipa±, 


(8.7) 


for ip G (ip±, V’i)- The complete action is S — Sf rcc + S int . 

The boundary equations of motion which follow from varying the action with respect 
to the fermion fields are, at r = 0: 


Vm- — e la i>~ = Y s ~ a+1 

ip + -e~"’lp- = yO +S+, 


(8.8) 


and the hermitian conjugate of these. The time derivative of the operators on the RHS of 
( |S.8| ) are determined by their commutation relations with the hamiltonian, O t O — —i[0 , H\. 
Using ( |8.7| ) and {S'-)., ST} = 1, one obtains 

d t (S-a+) = ^ (V>+ + a\ ^_) 

2 a _ ( 8 - 9 ) 
dt(a+S +) = (ip + + a\ ij)J) 
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and their hermitian conjugates (at r = 0). Combining (| 8 ]^), ( 8 T) and their hermitian 
conjugates, one obtains two independent equations at r = 0 : 


d t (V> 
dt (V> 


A 


+ - e 


^-) = — 7 - (V>+ + e ia ^_) 


A 2 


( 8 . 10 ) 


+ - e "if)-) = — (i/) + + e ia ^-) . 


The reflection S-matrices can be easily computed from (|8.10| ). The fermion fields have 
the mode expansions 

fM r °° 


= J d0 e~ B ! 2 ( 9)e~ ik ' r - Al^ (0)e ifc ' r ) 

_ rT/F r°o , x 

V>+M) - y— d0 e +e / 2 [A_(9)e- ikr + Al(9)e ik r ) , 


( 8 . 11 ) 


with ij)_ — ifj_ — i/j , . The variable 9 is the rapidity: k ■ r — M cosh(#)f — M sinh(0)r, 
and the A, are fermion annihilation/creation operators satisfying {A±(Q), A±(Q')} — 
5(9-9’). 

The scattering of the particles off the boundary at r = 0 can be formulated by formally 
introducing a boundary operator B satisfying the algebraic relational 


4m b = 4(»)4(-0) b. 


( 8 . 12 ) 


Separating the integral over 9 in ( |8.11| ) into 9 < 0 and 0 > 0. and making a change of 
variables 0 —> —9 in in the former, one obtains 


—e 2 ( iM cosh 9 + A+(Q) — e ia e 2 (iM cosh 9 — ] A^_(9) 


B 


A 2 


e 2 (iM cosh 9 + A\(—9) + e ia e 2 (iM cosh 9 — ) Alt (— 9) 


B 


e 2 (iM cosh 9 — A\_(9) — e ia e 2 (iM cosh 9 + ] A[_ (9) 


^ A 


—e 2 ( iM cosh 9 — J A J + (—9) + e ia e 2 ^iM cosh 9 + J All ( 


B 

A 2 A 


- 0 ) 


From this one can read off the reflection S-matrices: 


R+( 0 ) = 


sinh 9 


Rt(9) = 


R+(9 ) = RZ(9 ) = 


cosh((9 — 7 (9)) ’ 

cosh y( 0 ) 


sinh 9 


cosh(d — 7 (9)) 


B. 

(8.13) 


(8.14) 


cosh(d — 7 ( 0 )) ’ 


20 
















where 


(8.15) 


=7(0 — 


cosh# — iA 
cosh# + iA' 


A = 


A 2 


2 TO 


gap 


The parameter a is unphysical since it simply corresponds to a phase, hence we set it to 
zero. 

One can easily check that R satisfies the unitarity and crossing symmetry con¬ 
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R b J0)R£(-0) = Si 


Ri 


m 


2~ e =-*f T + ® 


in 


(8.16) 


where a — —a. 

One can also conhrm that these reflection matrices have the proper massless limit. 
For the left-movers, the massless limit is obtained by letting # —> # — a, and taking a —» oo, 
M —> 0 while keeping Me a /2 — /i held hxed. In this way, the massless dispersion relation 
([4.8|) is recovered from the massive one ([7.5|j. Taking this limit in R one hnds 


Rl 


R± 


# + #B 


with 


± tanh 

0, 


0 b 

ujb — fie — 


in 

T 


(8.17) 


A 2 


(8.18) 


This agrees with Q4. 1 1|) up to unphysical phases. The relation ( ^.18|) can now be used to 
express everything in terms of physical parameters. In particular, ( |8.18| ) combined with 

(0) £ ives 


A = 


UJq 


m 


(8.19) 


gap 


In terms of the physical energy 17, 


R+{E) = 


Rl = 


E 2 - 


m 


gap 


/4 


E 2 + D 2 /4 


E 


m 


gap 


E 2 ~ S?/4 + iDo v /B 2 -m | lp /4 
E 2 - wlH 


( 8 . 20 ) 


2 E 


E 2 - w 2 /4 + iuj 0 \JE 2 - m 2 ap /4^ 
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9. The Binding of Light to Atoms 


Using results from the last section, we now study the binding of light to atoms in 
photon bandgap materials in the simplest case of the free fermion point. Recall the physical 
parameters are the renormalized two-level splitting u3o, and the bandgap of the medium 

^gap- 


Light localized at, or bound to, the atom corresponds to a ‘boundary bound state’ of 
the kind described in |23[]. Boundary bound states exist if there are poles in the reflection 
S-matrices R b a {9 ) in ( fS.14|) on the physical strip 0 < Im 9 < ir. Poles in R occur when 
cosh(d — 7 (0)) = 0. Letting 0 — iu , poles occur at solutions of the equation 


sin 2 u + 2 A sin u + A 2 — 1 = 0. 


(9.1) 


Requiring the pole to be on the physical strip leads to a single boundary bound state at 
u — Ub satisfying 

sin it;, = 1 — A > 0. (9-2) 

Thus there exists a light-atom bound state if A < 1, or 


LOq A ^gap- 


(9.3) 


The necessary condition (|9.3|) for the existence of a bound state is simple to under¬ 
stand: the atom normally decays spontaneously by emitting photons of energy u)q, however 
such a photon cannot propagate in the medium and thus becomes bound to the atom. 

The energy of the boundary bound state above the ground state energy is M cos Ub, 
and this represents the binding energy: 


-E'bind — 


( TTl gap^O ( y ^0 




2 m 


gap 


1/2 


(9.4) 


The bound state is stable if the binding energy is less than the mass of the lightest particle, 

-S'bind T AT, Or Cdo <1 ^lgap/2- 

The existence of this bound state can be used to design a filter for light near the 
binding energy. Recall that with no atomic impurities, light with energy E < m gap is 
forbidden to propagate in the medium. The bound state is at an energy Eb < 'm gap / v^2 
due to the inequality ( |9.3| ), thus it is within the gap. If there are many atomic impurities, 
one can imagine that light of energy near Eb can resonate with the bound state and in 
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principle still propagate by tunneling along the impurities. One indication of this is the 
enhanced reflection amplitude for energies near E},\ 


R-AE) 


(m 


gap 


- Uo) U 0 


E — -Ebind 


4 El 


(9.5) 


One can determine the reflection S-matrices for particle scattering off of the bound 
state itself[^3|. In this simple case they turn out to be identical to the reflection S-matrices 
R b a for the scattering off the ground state computed in section 4. 


10. Conclusions 

We have found an exact solution to the problem of eigenstates of the atom-field inter¬ 
action in a spherical cavity which shows a rich spectrum in comparison with models based 
on the rotating wave approximation. The spectrum has solitons and bound states, in spite 
of the fact that there is no medium, i.e. we are dealing with a single atom in vacuum. In 
a sense our work shows that the well-known solitons of classical non-linear optics have a 
counterpart for a single atom in vacuum. 

By introducing a bulk mass term into the Kondo model in a way that preserves 
integrability we have shown how one can study analytically the occurance and properties 
of photon-atom bound states in photonic crystals. 

It is interesting to explore the possibility of utilizing the full spectrum of eigenstates 
in novel device applications. The solitonic mode being in a sense the most fundamental, 
one can imagine a ‘soliton laser’ based on the existence of this mode. 
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